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One of the most remarkable results to emerge from heavy-ion collisions over the past two decades is 
the striking regularity shown by particle yields at all energies. This has led to several very successful 
proposals describing particle yields over a very wide range of beam energies, reaching from 1 A GeV 
up to 200 A GeV, using only one or two parameters. A systematic comparison of these proposals is 
presented here. The conditions of fixed energy per particle, baryon+anti-baryon density, normalized 
entropy density as well as percolation model are investigated. The results are compared with the 
most recent chemical freeze-out parameters obtained in the thermal-statistical analysis of particle 
yields. The sensitivity and dependence of the results on parameters is analyzed and discussed. It is 
shown that in the energy range above the top AGS energy, within present accuracies, all chemical 
freeze-out criteria give a fairly good description of the particle yields. However, the low energy 
heavy-ion data favor the constant energy per particle as a unified condition of chemical particle 
l' freeze-out. This condition also shows the weakest sensitivity on model assumptions and parameters. 
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*^ 1 I INTRODUCTION 

H , 

One of the most remarkable results to emerge from relativistic heavy-ion collisions over the past years is the 
striking regularity shown by particle yields at all beam energies. From the lowest SIS up to the highest RHIC 
energies, all results on particle multiplicities are consistent with the assumption of chemical equilibrium in the 
final-state fireball produced after heavy- ion impact p|. The particle yields are found to be described, with 
remarkable precisio m by a thermal-statistical model that assumes approximate chemical equilibrium Q, 0, 0, 0, 
El H> Mi fiot ITU 1121 Il3l IT3. ITEj l . For a given collision energy, the thermal-statistical model with only two 
parameters, the temperature (T) and baryon chemical potential (/is), provides a very systematic description of 
particle yields. 

The partition function used in heavy-ion collisions has also been shown to be consistent with results obtained 
from lattice gauge theory (LGT) in the hadronic phase [l6l Il7). This, together with the phenomenological 
success of the thermal-statistical model, indicates that the resonance gas description of the hadronic phase is a 
remarkably accurate approximation of QCD thermodynamics in the hadronic confined phase. 

Using this description as a starting point, a further interesting systematic behavior of thermal parameters 
has emerged from particle yields in heavy-ion collisions from SIS up to RHIC. With increasing collision energy, 
there is an increase of the chemical freeze-out temperature, T, and a corresponding decrease of the baryon 
chemical potential, (1b ■ In the (T, (1b) ~ plane the freeze-out parameters lie on a curve connecting the lowest SIS 
through AGS, SPS and RHIC points, with a temperature at /is = that corresponds to the critical temperature 
expected for deconfinement in LGT 0, 0, ^ . Such a regular behavior of the freeze-out parameters has called 
for an interpretation based on the interaction dynamics. There are properties of the thermal fireball at chemical 
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freeze-out that are common to all collision energies. Such common properties provide unified chemical freeze-out 
conditions in heavy- ion collisions at all energies. 

The knowledge of such unified freeze-out properties is useful since: 

• they make it possible to determine particle excitation functions and to study their systematic properties 
|20j . Consequently, general conclusions can be drawn about the energy dependence of particle production 
in heavy-ion collisions; 

• the freeze-out conditions are relevant for a successful description of particle spectra; 

• on the theoretical level the unified freeze-out conditions are relevant for the understanding of particle 
productio n dy namics. The chemical freeze-out line has been interpreted as the inelasticity line in heavy- ion 
collisions plf . as it separates the stage where the evolution dynamics of the collision fireball is dominated 
by inelastic instead of elastic processes. 

The first unified chemical freeze-out criterion was recognized |22ll23| when comparing the thermal parameters 
at SIS energy with those obtained at SPS. It was shown that the average energy per particle at SIS energy reaches 
approximately the same value of 1 GcV as calculated at the critical temperature expected for deconfinement at 
fiB = 0. In addition, known results for chemical freeze-out parameters at the AGS also reproduced the same 
value of energy per particle. Thus, it was suggested that the condition of a fixed energy per hadron is the 
chemical freeze-out criterion in heavy-ion collisions. More precise determinations prefer a value (E) / (N) sa 
1.08 GeV. 

In addition to the fixed (E) / (N) criterion, alternative proposals have been made to describe chemical freeze- 
out in heavy- ion collisions at all energies pi El Ej I27L I2I El l30| : 

• a fixed value for the sum of baryon and anti-baryon densities, ub + of approximately 0.12/fm 3 |24j |: 

• a self-consistent equation for the densities based on geometric estimates using percolation theory p5| : 

• a fixed value of the entropy density, s/T 3 , of approximately 7 prl I27H2M I29L l30| . 

All of these proposals have been used with considerable success, and it is the purpose of the present article 
to study and compare these different criteria for particle freeze-out to constrain their validity on the basis of 
available experimental data and corresponding statistical model analyses. It is not always straightforward to 
judge their quality, since small variations in parameters can sometimes lead to substantial shifts in results. In 
our analysis, we quantify the explicit dependence of different freeze-out criteria on excluded volume corrections, 
the resonance composition of the fireball, strangeness content and the corrections related to the canonical 
implementation of strangeness conservation. In particular, excluded volume corrections |3lL l32| have a very 
strong impact on the validity of most of the criteria. 

We show that in the energy range above the top AGS beam energy and within present accuracies, it is rather 
difficult to distinguish between different freeze-out conditions. However, the lower energy heavy-ion data favor 
the constant energy per particle as the unified condition of chemical particle freeze-out. It also exhibits the 
least sensitivity to model assumptions and the most stability to small changes in the values of the parameters. 

II. GLOBAL DESCRIPTION 

The results from particle yields obtained by different groups for the thermal parameters are not always 
identical. The differences are small and never exceed the level of a few percent. Nevertheless, such differences 
play a role when one tries to determine unified descriptions over a wide range of beam energies. The 
origin of these differences has not been fully studied until now. However, to a large extent these differences 
could be attributed to different assumptions on the number of resonances, their widths, the treatment of 
weak decays, the description of strangeness saturation and how, or if at all, repulsive interactions have 
been included in the calculations. To pin down the differences would involve a very detailed knowledge of 
programming codes. This is not yet possible, as codes are not always made available. In order to have 
a consistent presentation, we have made some choices which are biased by our own experience: the code 
used in this paper is available for inspection [33j. Our analysis relies on the most recent results obtained in 
statistical-thermal model fits to Au+Au and Pb+Pb systems, performed by numerous groups over a wide 
range of energies 0, IE Hi U EH El El EL El ■ These results are summarized in Table Q] and are also 
included in Figure ^ 
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Collision System and Energy 


Ref. 


T 


/US 


Include in 






(MeV) 


(MeV) 


Fits 


RHIC 


Au+Au y/s~ = 200 AGeV 


m 


177+7 


29+6 






M 


163+4 


24+4 






[15J 


165.6+4.5 


28.5+3.7 




Au+Au y/s = 130 AGeV 


® 


169+6 


38.1+4.2 


/ 




® 


174+7 


46+5 


/ 




[14] 


165+7 


41+5 


/ 


SPS 


Pb+Pb 158AGeV f = 17.3 AGeV 


m 


157.5+2.2 


248.9+8.2 


/ 




m 


154.6+2.7 


245.9+10.0 


/ 






161.0+6.0 


260.0+30 


/ 


Pb+Pb 80AGeV yfs = 12.3 AGeV 


® 


153.5+4.1 


298.2+9.6 


/ 




m 


149.9+5.1 


293.8+11.0 


/ 




[12] 


155.0+5.0 


284.0+15.0 


/ 


Pb+Pb 40AGeV yfs = 8.77 AGeV 


m 


146.1+3.0 


382.4+9.1 


/ 




w 


143.0+3.1 


380.8+8.9 


/ 




[12] 


148.0+5.0 


367.0+14.0 


/ 


Pb+Pb 30AGeV yfs~ = 7.62 AGeV 


w 


140.1+3.3 


413.7+16.3 


/ 




m 


144.3+4.7 


406.0+19.1 


/ 


Pb+Pb 20AGeV yfs = 6.27 AGeV 


m 


131.3+4.5 


466.7+12.9 


/ 




m 


135.8+5.2 


472.5+13.7 


/ 


AGS 


Au+Au 11.6AGeV yfs = 4.86 AGeV 


m 


118.7+3.1 


554.4+13.0 


/ 




w 


119.2+5.3 


578.8+15.4 


/ 




[12J 


123.0+5.0 


558.0+15.0 


/ 


SIS 


Au+Au l.OAGeV yfs = 2.32 AGeV 


m 


52+1.5 


822 


/ 




m 


49.6+1 


810+15 


/ 






49.7+1.1 


818+15 


/ 






58+4 


792+7 


/ 


Au+Au 0.8AGeV yfs = 2.24 AGeV 


[101 


54+2 


808+5 


/ 



TABLE I: Results obtained in statistical-thermal model fits to Au+Au and Pb+Pb collision systems by numerous groups 
over a wide range of energies. The checked entries have been included in the fits to determine the T(/j,b) and ^s( v / s) 
parameterizations (Equations Q and respectively). 



Before discussing regularities in Table [I] we present first a polynomial fit using the checked entries in Table |U 
as input. This gives, 

T(fi B )=a-bfj, 2 B -cfi%. (1) 

where a = 0.166 ± 0.002 GeV, b = 0.139 ± 0.016 GeV^ 1 and c = 0.053 ± 0.021 GeV~ 3 . The yfs = 200 AGeV 
points were not included in the determination of Eq. Q). The comparison of the above equation with the 
•^/s-dependent chemical potential from Table [I] is shown in Figure ^ It is also straightforward to verify that at 
yfs = 200 GeV Eq. (JTJ is consistent within errors with the expected results from thermal model analyses. The 
energy dependence of the baryon chemical potential can be parameterized as |34j : 

^ = rrbr (2) 

with d = 1.308±0.028 GeV and e = 0.273±0.008 GeV" 1 . The temperature T and the baryon chemical potential 
Hb are shown in Figure [3 as a function of the beam energy. It provides a good quantitative description of the 
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H R (GeV) 



FIG. 1: Values of \ib and T for different energies. The solid line is a parameterization corresponding to T(^ib) ~ 
0.17 - 0.13^1 - 0.06^1 (see text). 

corresponding thermal model results. Eqs. JU and J5J) lead to the energy dependence of the chemical freeze-out 
temperature and baryon chemical potential shown in Figure [3 In heavy ion collisions at the LHC Eq. (0 leads 
to a value /j,b ~ 1 MeV. 

The parameterization given above should of course be used with caution, as small changes in the parameters 
lead to a significant change in the energy dependence, since particle yields depend on /j,b and T exponentially. 
A particularly 

Having established a basic parameterization of freeze-out conditions at different energies, we now proceed to 
compare various unified conditions for the chemical freeze-out curve. In the next section we will discuss how 
the freeze-out criteria mentioned in the introduction live up to a comparison with the description presented 
here. We will, in particular, discuss the dependence on parameters such as the strangeness suppression factor 
7s, canonical corrections to strangeness conservation, the cut-off of the hadronic mass spectrum and, finally, 
excluded volume effects. Clearly, a superficial examination shows that many criteria can be remarkably effective 
at describing chemical freeze-out in the restricted energy range. The largest discrepancies, as it will be argued, 
appear at low energies. However, a closer examination reveals that some criteria are more robust than others 
when small changes in parameters are allowed. 

III. DESCRIPTION OF FREEZE-OUT CRITERIA 

All criteria considered here are based on the hadron resonance gas model that was successfully used to describe 
particle yields and their ratios in heavy-ion collisions. Although this model is well known in the literature, for 
completeness we will summarize the basic ideas and formulae required to quantify each freeze-out condition. 

The model is defined by the partition function of a resonance gas with free particle dispersion relations for 
all constituents The resulting observables, however, account for interactions between hadrons, to the extent 
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that the thermodynamics of an interacting system of elementary hadrons can effectively be approximated by a 
mixture of ideal gases of stable particles and resonances. Thus, the thermodynamical observables are expressed 
as a sum over all mesonic and baryonic degrees of freedom, as well as their known resonances. 

For our further discussion, the relevant quantities are: the total particle density n, total density of baryons 
iib + Tig, as well as energy e and net baryon number ns — ng densities. All of these observables are functions 
of temperature T and chemical potentials /is,/XQ and fis related to conservation of baryon number, electric 
charge and strangeness, respectively. The last two potentials are fixed from the initial conditions. Thus, we are 
left with only two independent variables: T and [Lb- The results for all relevant quantities obtained from the 
hadron resonance gas partition function are as follows: 
The particle density is given by, 

n = £»/^r/«. (3) 

the entropy density by, 

3 = E A / P" d P ["/< ln fi T (1 T fi) ln(l T /*)] , (4) 
< 27r Jo 




FIG. 2: Energy dependence of the chemical freeze-out parameters T and fj,B- The curves have been obtained using a 
parameterization discussed in the text. 
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and, finally, the energy density, 
where, 

with the upper (lower) signs being for Fermi-Dirac (Bose-Einstein) statistics. The total density of baryons is 
calculated as in but restricting the sum to baryons ng and their antiparticles rig. All of these results 
correspond to the grand-canonical formulation of the conservation laws. For central heavy-ion collisions at 
energies below AGS, strangeness conservation on average is not adequate and has to be implemented exactly in 
terms of the canonical ensemble || . A more general result for basic thermodynamic quantities in the canonical 
formulation of strangeness conservation can be found in Ref. [p. T he thermal particle phase-space of hadronic 
resonances has also to be modified by including finite width yj, l35| . This modification amounts to convoluting 
the particle momentum distribution in Eqs. J3][5j) with the Breit-Wigner resonance form-factor. 

The strange particle phase-space in Eqs. is controlled by the 7 S parameter j3(J. The value of y s = 

1 corresponds to a chemical equilibrium distribution, whereas 7 S ^ 1 indicates a deviation from chemical 
equilibrium. In addition, Eqs. are only valid for point-like particles. Considering hadrons as extended 

objects would require an implementation of interactions. This can be done in the hard-core |32j| and self- 
consistent implementation [3l]|. The first approach amounts to the multiplication of all results in Eqs. ©, (0} 
and JSJ by a factor, 



F(T,hb,V )= -, (7) 

1 + Von{T,^ B ) 

where the repulsive strength is determined by the volume parameter Vq. 

In the self-consistent implementation of the repulsive interaction there is, in addition, a shift in chemical 
potential that guarantees that all observables Eqs. ©-© are obtained as corresponding derivatives from the 
partition function |3lT ] . The explicit results and a detailed description can be found in Ref. |3l| . 

Before applying the above results to quantify different freeze-out criteria and comparing their predications 
with the thermal model analysis of heavy-ion data, we first briefly summarize basic concepts of these criteria. 



A. Percolation Model 

Percolation theory was successfully applied to the description of critical properties of QCD matter 133 . 
In particular, color deconfinement in pure gauge theory can be treated as a percolation phenomenon [38j. 
The models based on percolation theory also provide a qualitative description of charmonium production and 
suppression in heavy-ion collisions |3^ |. 

In Ref. |25| , percolation theory was used to formulate and quantify chemical freeze-out conditions in heavy- ion 
collisions. It is assumed that, in the sector of vanishing baryon density, hadronic matter freezes out according 
to the resonance gas approximation and by vacuum percolation, while the region of finite baryon density freezes 
out according to ba ryon percolation. The condition that describes the freeze-out line in heavy-ion collisions has 
been formulated as 25] : 



n(T, n) 



1.24 



n(T,n) 



0.34 



nB(T,n) 
n(T,n) 



(8) 



The volume Vu is the hadronic size, which corresponds to a radius of approximately w 0.8 fm. The numbers 
1.24 and 0.34 appearing in Eq. JHJ are calculated within percolation theory ^4|. The smallest one is determined 
when the size of the largest cluster falls below the size of the overall spatial volume. On the other hand, the 
number 1.24 is determined by the disappearance of any large-scale vacuum, and only the strongly interacting 
medium spans the entire space. 

Eq. JSJl , together with the results of the statistical-thermal model of the hadron resonance gas on the total 
density of all hadrons n(T, iib) and baryon number density ns(T, //a), provide a unique solution T = T(fis) 
that characterizes the chemical freeze-out line in the {T^b) ~ plane within the percolation model. 



7 



B. Fixed Baryon+Anti-Baryon Density 

The contributions of baryons to thermodynamics is usually characterized by the net baryon density, i.e. the 
difference between the density of baryons ng and anti-baryons ng. In heavy-ion collisions, this varies strongly 
with collision energy. At SIS it is ~ 1/3 of nuclear matter density (po — 0.17/fm 3 ) and it strongly increases to 
reach a maximum at top AGS energies. It then gently decreases up to RHIC, where it is only ~ 1/10 of po- 
Due to baryon number conservation (B — A part ), the variation of ng — ng = B/V reflects the corresponding 
change in the volume of the fireball. Indeed, a comparison of the total yields of different particle species in 
heavy-ion collisions with calculations in the statistical model, or, alternatively, u sing two pion correlations, has 
shown that the volume has a minimum between AGS and the lowest SPS energy ;40j and exhibits a dependence 
on \/s opposite to that of ng. 

Despite this large variation, it was first noticed in Ref. p4| that the sum of baryon and anti-baryon densities, 
n B + n B , remains remarkably constant. Thus, it was pr~osed that in heavy -ion collisions the chemical 
freeze-out curve can be obtained from the condition, 

n B + n B ~ 0.12fm~ 3 , (9) 

of fixed density of total number of baryons and anti-baryons. 

The observation |40g that Eq. JSJ) is satisfied in high energy heavy-ion collisions was important phenomeno- 
logically to explain the equivalence at 160 and 40 AGeV of the enhancement of the low mass dilepton yield 
observed by the CERES Collaboration in Au-Pb collisions j^jj and recently by NA60 in In-In reactions 
Theoretically, the enhancement was argued ^3| to be due to rho meson broadening through its rescattering 
with surrounding baryons. Thus, the broadening is determined by the total density of baryons. The fact that in 
Au-Pb collisions at 160 and 40 AGeV the total baryon density is comparable, gave an argument for the observed 
similar structure of the low mass dilepton yield at these energies. 

Clearly, in low energy heavy-ion collisions between SIS and AGS, the production of anti-baryons is suppressed. 
Thus, in this energy range the net baryon number density, to a good approximation, is equivalent to the total 
density of baryons, that is ng + rig . Consequently, since the volume at freeze-out from SIS to AGS extracted 
from HBT correlations 40] drops by a factor 1/3, there should be a corresponding increase in ng +n§. Because 
of the above, the freeze-out condition © is not expected to work in low energy heavy-ion collisions. The 
question remains as to at what energy the violation of 10 as the condition of particle freeze-out sets in. The 
answer to this question requires a quantitative analysis, which is given in the next section. 



C. Fixed Entropy Density over T 3 

The chemical freeze-out condition of fixed entropy density over temperature to the third power, 

^3- 7 > (10) 

was proposed independently from two very different considerations. In Ref. pfll l30j | it was used to extrapolate 
lattice gauge results from pg = to finite values of pg by keeping s/T 3 fixed. In Ref. |2fll27j the entropy density 
was used to separate between a baryon-dominant region and a meson-dominant one, in order to understand the 
rapid change in certain particle ratios observed at lower SPS energies by the NA49 collaboration [45|. In the 
course of this analysis, it was noted that, over a broad energy range, the total entropy density divided by T 3 
remained constant, despite very large changes in the baryon chemical potential and the temperature. 

The freeze-out condition l|10[l is rather surprising. In the ideal gas of massless constituents and for a net 
baryon free system, the s/T 3 ratio describes the number of degrees of freedom. In the case of massive particles 
and at finite chemical potential, this ratio can be parameterized as s/T 3 ~ a(T,pg). Thus, for fixed T and 
pg, the parameter a(T,pg) describes the effective number of degrees of freedom in a system. The condition 
(|10|l tells us that this effective number of degrees of freedom a(T,pg) is common for all collision energies. In 
the thermodynamical context, this also means that a(T,pg) is independent of thermal parameters at chemical 
freeze-out. This is quite unexpected, as changing the collision energy implies not only changing the thermal 
parameters but also the hadronic composition of the fireball. Naively one would expect that in such a case the 
number of degrees of freedom a(T,pg) should be a strongly varying function of collision energy. 
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D. Fixed (E) / (N) ~ 1 GeV 

This criterion was first proposed in Ref. [22l l23j as a link unifyingresults obtained at SIS with top AGS and 
SPS. Since then, the description has proven its validity also at RHIC ]E@ an d at the lower SPS energies Q,0,E3- 
This condition was successfully used in the literature to make predictions |46| of freeze-out parameters at SPS 
energies of 40 and 80 AGeV for Pb-Pb collisions long before the data were taken. These predictions turned out 
to be in agreement with the statistical model analysis of recent experimental data obtained at these collision 
energies Q- 

The reason why (E) / (N) is constant in heavy-ion collisions from SIS up to RHIC can be intuitively explained 
as follows: In a non rclativistic system (E) / (N) ~ (to) + |T, with (m) being the average mass in the thermal 
fireball |47|. At SIS, the thermal mass is of the order of the nucleon mass, whereas the freeze-out temperature 
T f ~ 50 MeV. Thus, at SIS the (E) / (N) ratio is indeed comparable with 1 GeV. At SPS and RHIC, the 
leading particles in the final state are pions. However, at chemical freeze-out most of the pions are still hidden 
in the mesonic and baryonic resonances. Thus, here the average thermal mass corresponds approximately to 
the /9-meson mass. Consequently, with Tf ~ 160-170 MeV a comparable value of (E)/(N) is obtained at SPS 
and RHIC as in a much lower energy at SIS. 

We have already indicated that the chemical frcczc-out line dynamically coincides with the inelasticity line in 
heavy-ion collisions. Thus, also the (E) / (N) ~ 1 GeV condition should have the above dynamical interpretation. 
Recently, the correlation between the (E) / (N) ~ 1 GeV condition and inelasticity has been investigated in 
central Pb-Pb collisions at the SPS in terms of the Ultra-relativistic Quantum Molecular Dynamics model 
(UrQMD) A detailed study has shown that there is a clear correlation between the chemical break-up in 
terms of inelastic scattering rates and the rapid decrease in energy per particle. If (E) / (N) approaches the value 
of 1 GeV the inelastic scattering rates drop substantially and further evolution is due to elastic collisions that 
preserve the chemical composition of the collision fireball. Following the above UrQMD results, it is conceivable 
that (E) / (N) ~ 1 GeV is related with inelasticity in heavy- ion collisions. This could dynamically justify the 
applicability of this freeze-out criterion. 

It was noted recently in Ref. that the results from RHIC favor a higher value for (E) / (N). This 
was interpreted as a change in the baryonic composition of a two-component thermal source. We find that a 
consistent description using (E) / (N) = 1.08 GeV is quite adequate and there is no need for different values at 
different energies. The comparison of each of these freeze-out criteria with all known freeze-out parameters is 
shown in Figure 03 

IV. SENSITIVITY AND ROBUSTNESS OF FREEZE-OUT CRITERIA 

The chemical freeze-out criteria discussed in the last section are not quite model independent. Their pre- 
dictions are sensitive to various assumptions and to the choice of parameters. In this section, we present a 
detailed discussion of the robustness of these criteria. In this context, we consider the influence of the repul- 
sive hadronic interactions, strangeness suppression and hadron mass spectrum cut-off on the predictions of the 
different frcczc-out criteria. 



A. Sensitivity on excluded volume corrections 

The (E) I (N) criterion involves the ratio of two extensive quantities and, hence, is nearly independent of 
excluded volume corrections. In the Boltzmann approximation, these corrections cancel exactly, provided that 
the same proper volume parameter is applied. Otherwise, for example, pions and protons pick up different 
correction factors through the shift in chemical potentials. In quantum statistics, the shift in chemical potentials 
always leads to a (small) change. The energy per particle is also not strongly influenced if repulsive interactions 
are implemented through a mean-field potential. The freeze-out condition of fixed entropy density divided by 
T 3 , as well as that of fixed total density of baryons plus anti-baryons, use extensive thermodynamic variables. 
Hence, the dependence on the excluded volume corrections is, in general, not negligible. However, as it turns 
out, only a moderate dependence is found for the constant ng + rig criterion, which is shown in Figure^] The 
repulsive strength of hadrons is parameterized there by the radius parameter. The freeze-out curve moves to 
higher values of T as the radius increases. A change in radius of 0.1 fm easily leads to a change in temperature 
which can be larger than 10 MeV. The repulsive interactions suppress thermal particle phase-space. Thus, to 
reproduce the same value of total density of baryons, one needs to increase temperature at fixed hb- 
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FIG. 3: Description of chemical freeze-out by four different criteria discussed in the text. 

The percolation model criterion is determined by the total particle and the net baryon number density. It thus 
shows a similar dependence on excluded volume radius as discussed above. This model is based on geometric 
considerations which are affected by volume parameters. The changes are almost identical to the previous case 
and are not shown explicitly. 

Contrary to the percolation and fixed baryon density conditions, the constant s/T 3 criterion shows a wild 
dependence on the excluded volume radius (as seen in Figure [Q. This casts considerable doubt on the value of 
this criterion for model building. It is clear that only the complete neglect of excluded volume corrections leads 
to a realistic freeze-out curve. It is well known, however, that such corrections cannot be ignored, especially 
at low temperatures, as this leads e.g. to a totally unrealistic picture for the phase transition and is also 
inconsistent with the properties of the nuclear potential. 

B. Sensitivity on resonance spectrum mass cut-off 

The thermal-statistical model relies very strongly on input from the particle data table • Heavy resonances 
are, in general, suppressed by the Boltzmann factor e~ E l T . However, when the temperature or baryon chemical 
potential increase, the contribution of heavy resonances also increases and this becomes especially noteworthy 
in the small \xb - high T region. As the information about heavy resonances is fragmentary, several analyses 
have introduced a cut-off in masses, such that resonances heavier than this cut-off are not taken into account 
in the thermal-statistical model. The effects of this mass cut-off are similar for all criteria considered. As an 
example, we show the influence on the (E) / (N) — 1 GeV criterion and on the finite density ns+ng = 0.12/fm 3 
criterion in Figure [3] It is seen in Figure [S] that the largest influence of the mass cut in the hadron spectrum 
appears at high temperature and low fis- At fixed hb, the changes in the freeze-out temperature never exceed 
10 MeV. 

The cut-off on the hadronic mass spectrum in the thermal-statistical model is particularly relevant when 
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FIG. 4: The effect of excluded volume corrections on the constant ub + n B (bottom) and constant s/T 3 (top) freeze-out 
criterion. 

calculating particle yields. This is because the decay properties of heavier resonances into lighter hadrons are 
needed and are not always well established. However, for the description of basic thermodynamic observables 
characterizing the properties of the hadronic medium, like the energy, entropy, particle and baryon number 
density, the only resonance characteristics needed are the masses, degeneracy factors and quantum numbers. In 
most cases, this information is known even for heavy mesonic and baryonic states. 



C. Sensitivity on strangeness suppression 

Strange particle yields, their ratios and momentum distributions are relevant observables that characterize 
the origin, composition and production dynamics in heavy- ion collisions |ll l36| . However, their contribution 
to global thermodynamic observables, summarized in Eqs. I|3l5f) . is not essential. In low energy heavy- ion 
collisions at SIS up to AGS, strange particles are only rarely produced and can be almost completely neglected 
when considering global thermodynamic characteristics of the collision fireball. At SPS and RHIC, due to the 
much higher temperature, the strangeness degrees of freedom contribute about 25% to the total particle density. 
Due to this, one should not expect very large changes in the predictions of different freeze-out criteria when 
modifying the strange particle content. 

In the statistical model under consideration, a modification of the strange particles from chemical equilibrium 
is described by the 7 S factor or by the canonical constraints imposed by exact strangeness conservation. 
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FIG. 5: The effect of the resonance spectrum mass cut-off on the chemical freeze-out curve: top constant (E) / (N), 
bottom ub + ng freeze-out criterion. 



The strangeness-suppression (or enhancement) factor 7 S modifies the multiplicities of strange and anti-strange 
hadrons, parameterizing a deviation from chemical equilibrium, in the Boltzmann approximation, in the manner 
given by the equation below, 

(«*) = lT ("^equilibrium ■ ( U ) 

It does not affect the non-strange hadrons except for those contributions coming from decays of strange reso- 
nances. Its influence on the chemical freeze-out curve is therefore small and, at low energies, negligible. This 
is shown for the s/T 3 criterion and for the percolation model in Figure |SJ The effect on the other freeze-out 
criteria is similar and not shown explicitly. 

For small systems, various corrections become important. In our considerations, these affect mainly the 
production of strange particles in heavy-ion collisions at low energies due to the very small number of produced 
particles. In this case, strangeness conservation should be treated exactly, which results in so-called canonical 
corrections. These depend on the volume over which strangeness is balanced. This correlation volume does not 
need to be as large as the volume of the overall system. Several calculations have been made where strangeness 
conservation is imposed over volumes which are smaller than the overall volume. These are parameterized by a 
canonical radius. In Figure one typical example shows the influence of this canonical radius on the freeze-out 
criteria. Similarly as for the strangeness suppression factor, since the relative number of strange particles is 
small, the effect of the canonical radius on the chemical freeze-out curve is small. All criteria are only weakly 
dependent on variations in the canonical parameter. Small differences are only seen above AGS energies, where 
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the contribution of the strangeness sector is not negligible. At low energy, this effect is entirely irrelevant. 
However, for the discussion of strange particles or their yields, the influence of canonical effects is essential 
|8l llll lift Ifilj . This is particularly the case for multistrange baryons in low energy central collisions, as well as 
in peripheral A-A collisions at higher energies [i^. l5lj . The 7 S factor has a similar influence on strange and 
multistrange particle yields 0, Q • 

Comparing different freeze-out criteria, we have been bound by the knowledge of the lowest chemical freeze- 
out parameters at SIS energy. Clearly, the thermal-statistical model considered here is limited and cannot 
be used if the assumption of energy conservation on the average is not valid. This indeed is the case in low 
energy heavy-ion collisions, where the full micro-canonical formulation |52| of the partition function is required. 
This also means that the extrapolation of different freeze-out criteria to the T ~ limit and within the grand- 
canonical formulation is questionable. Until now the model has shown its validity in the description of heavy-ion 
data up to beam energies of about 800 MeV. 

There is a class of further criteria (e.g. the constant energy or total particle number density) that are 
consistent with the freeze-out parameters in the energy range above the lowest SPS energies. However, since 
these conditions are based on intensive observables, they suffer from variations with the parameterization of 
repulsive interactions. In addition, both of these criteria are inconsistent with the SIS freeze-out parameters. 
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FIG. 7: The effect of the canonical radius on the constant s/T 3 freeze-out criterion. 

V. SUMMARY AND CONCLUSIONS 

From the results presented above one can conclude that, on a phenomenological level, the various freeze- 
out criteria give quite satisfactory descriptions of the particle multiplicities measured in heavy-ion collisions. 
Deviations can be reduced by suitably modifying input parameters in the different models. However, it is clear 
that some of the criteria are not very stable or robust when one applies small changes to the input considerations 
and are, therefore, more model dependent. 

A further method to differentiate freeze-out conditions would be to calculate particle excitation functions. 
The differences between such criteria would be even more transparent, as some particle excitation functions are 
very sensitive to the value of the thermal parameters. A good example is the K/tt ratio which at SIS energy is 
a very sensitive probe of temperature and canonical suppression effects. The energy per particle criterion has 
been shown to provide a consistent description of the particle excitation functions over the whole energy range. 
Some of the results of this criterion, like the energy dependence of the A/ir and K + /K~ ratios for different 
centrality and collision energies, have been confirmed by recent heavy-ion data. 

The main differences between the chemical freeze-out criteria analyzed here appear in the region of low 
temperature and large baryon chemical potential, that is for SIS up to AGS energy. This region obviously 
favors the (E) / (N) criterion as the most appropriate one to describe chemical freeze-out over all energies. This 
is supported quantitatively by the fact that this freeze-out condition has the lowest x 2 P er degree of freedom 
of all the criteria considered in this analysis and was shown to provide a good description of different particle 
excitation functions in heavy-ion collisions. The energy per particle is also consistent through transport model 
calculations with the dynamical interpretation of chemical freeze-out as the inelasticity condition in heavy-ion 
collisions. 
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